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The noncentral chi-square distribution occurs in noise interjerence 
problems. When the number of degrees of freedom becomes large, the middle 
portion of the distribution is given by the central limit theorem, and the tails 
by a classical saddle point expansion. Here recent work by N. Bleistein and 
F. Ursell on "uniform" asymptotic expansions is combined and extended 
to obtain an asymptotic series which apparently holds over the entire range 
of the distribution. General methods for expanding saddle point integrals in 
uniform asymptotic series are discussed. Recurrence relations are given for 
the coefficients in two typical cases, (i) when there are two saddle points and 
(ii) when there is only one saddle point but it lies near a pole or a branch 
point. 

I. INTRODUCTION 

This paper deals with the problem of obtaining asymptotic scries 
for the complex integral 

J = f t^g(t) exp [xh(t)] dt (1) 

when x becomes large. Problems of this sort are quite often encoun- 
tered in applied mathematics, particularly in wave propagation. The 
material presented here grew out of some recent work by G. H. Rob- 
ertson 1 on the "Marcum Q-Function." This function, which appears 
in the study of radar interference, gives the distribution of the random 
variable (noncentral x~) 

z = (l/x) ± yl . (2) 

n-l 

1971 
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Here £ is a positive integer and j/i , y$ , . . . , y w are independent gaus- 
sian random variables with unit variances and mean values which 
may be different. 

Mr. Robertson has devised an algorithm for computing the Q-func- 
tion which may be used for a wide range of the parameters appearing 
in the function (that is, in the noncentral x 2 distribution). In an 
earlier paper on information theory, by working with an integral of 
the type in equation (1), I had obtained an asymptotic (for large 
x) expression for the tails of the distribution. 2 However, comparison 
with results obtained by Robertson showed that my expression failed 
badly in the central part of the distribution where the central limit 
theorem holds. 

The need for an asymptotic expansion which holds uniformly over 
the entire range of the distribution led to a study of the recent work 
on "uniform" asymptotic expansions of integrals. The first part of 
this paper is an exposition, plus extensions and generalizations, of 
some of the procedures which have been used to obtain uniform 
asymptotic expansions of integrals of the type in equation (1). The 
theory is then applied to the noncentral x 2 distribution. 

Two procedures are considered. For convenience, we call them the 
"Bleistein method" 3 and the "Ursell method." 4 Although these names 
are among the best that suggest themselves, they are not entirely 
satisfactoiy because they contain no hint of the earlier work by 
others, especially Olver, Chester, Friedman, and Ursell. B - 8 Here we 
have recast the underlying ideas used by Bleistein and Ursell into 
forms better suited to our purpose. 

Both methods lead to the same asymptotic series. The Bleistein 
method gives a compact expression for the coefficients in the expan- 
sion. However, from the few examples that have been studied, it 
appears that the labor required to reduce this compact expression to 
a computable form is at least as great as that required by the Ursell 
method. 

Section III and Appendices A, B, and C are concerned with a prelimi- 
nary change of variable in the integral J. The case, denoted by "X =1" 
for brevity, in which the exponent X is a positive integer, is discussed in 
Sections IV, V and VI. This material is applied to the problem of two 
saddle points in Appendix E. The case in which X is general, denoted 
briefly by "X j* 1," is discussed in Sections VII, VIII, and IX, and in the 
examples in Appendices F, G, and H. The results of Section IX are ap- 
plied in section X to obtain the desired type of expansion for the non- 
central x distribution. Useful results regarding classical saddle point 
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expansions are stated in Appendix D. Some of the results given in 
Appendix F for the general case of a saddle point near a branch point are 
applied in Appendix G to obtain an asymptotic series for the Poisson- 
Charlier polynomial, a polynomial of interest in traffic theory. 

II. STATEMENT OF PROBLEM 

The general problem is to obtain an asymptotic series for the inte- 
gral J defined by equation (1) when x becomes large and most of the 
contribution to J arises from a (rather loosely defined) "critical re- 
gion" around t = 0. The path of integration U is supposed to start 
and end at \t\ = oo in "valleys" in the complex i-plane where 
exp[xh(t)] -> as \t\ -» oo. Let the starting and ending valleys be 
denoted by S and E, respectively. The path U starts in S, climbs up 
to and passes through the critical region, and then descends down 
into E. 

The functions h(t) and q{t) are analytic in the critical region; 
and one or more saddle points, that is, points where h'(t) = dh(t) /dt 
vanishes, lie in the critical region. We assume h(0) = and that x is 
real and positive. If x were complex, the factor exp(i arg x) could be 
included in h(t). 

The path IJ may be deformed into a path D consisting of (i) paths 
of steepest descent which pass through some or all of the saddle 
points plus possibly (u) loops around branch cuts and poles. The 
path D is independent of x. When x is extremely large, all but a 
negligible part of J arises from contributions of veiy small portions 
of D. If t = is a singularity, one portion may lie close to t = 0. 
Another portion is centered on the highest (that is, largest exp[xh(t) ] ) 
saddle point. If the two highest saddle points are of the same height, 
a portion is centered on each, and so on. Thus when x is extremely 
large, the asymptotic scries for J may be obtained by the classical 
or "usual" saddle point method. 

However, we may wish to compute J for values of x which, though 
large, are not large enough to allow J to be evaluated by the classical 
saddle point method. For such x's the highest saddle points and the 
singularity (for X ^ 1) at t = cannot be treated separately, that is, their 
interaction must be taken into account. If other saddle points of lesser 
height he in the critical region, they must also be considered. This is the 
range of x of interest here. Our problem is to obtain the appropriate 
expansion of J in descending powers of x. The type of expansion we seek 
is shown in equation (46) for J. 
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This type and the type shown in equation (15) have occurred in 
earlier publications 3-8 and have been called "uniform" asymptotic 
expansions because they hold uniformly as a saddle point approaches 
a singularity or another saddle point. 

III. CHANGE OF VARIABLE 

In Appendix A it is shown that, in the critical region, h(t) behaves 
much like a polynomial of degree p + 1 in t. Here p is the number of 
saddle points in the critical region. This suggests changing the variable 
of integration from t to v where 

F(v) = h(t) (3) 

and F(v) is a polynomial of degree y. + 1 in v. When F(v) is known, 
solving (3) for v as a function of t gives p + 1 branches. The branch 
chosen for the change of variable is the one for which dt/dv « c through- 
out the critical region, c being a constant. That one and only one of the 
n + 1 branches has this property is rendered plausible by the discussion 
in Appendix A. 

Fortunately we do not have to solve equation (3) to obtain the 
asymptotic series we desire. However, for some steps we do need the 
values of dt/dv and higher derivatives at the saddle points. These 
may be obtained by repeated differentiation of (3) . 

F(v) is not uniquely determined by h(t). The factors which influ- 
ence its choice are reviewed in Appendix B. 

The change of variable from t to v carries the integral (1) for J into 

J = J /"'/(y) exp [xF(v)] dv (4) 

where 

j{v) = g(t) {t/vf-H n \ t a) - dt/dv. (5) 

The path of integration L starts in the v-plane valley corresponding 
to valley S in the £-plane, passes through the critical region surround- 
ing v = 0, then descends into the f-plane valley corresponding to 
valley E. 

IV. THE BLEISTEIN METHOD FOR A = 1 

For the case A = 1, the integral J becomes 
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- f g(t) exp [xh(t)] dt = / /(«) exp [»?(»)] dv, 



(G) 



m = g(t)f v = g(¥ ,) 



The Bleistein method begins by constructing a polynomial P {v) of 
degree p - 1 such that P (v r ) = f(v r ), r = 1, 2, . . . , p where Vi, v 2 , 
. . . , Vp are the zeros, assumed simple, of F'(v). By Lagrange's inter- 
polation formula, 

P (v) _ y mm (7 ) 

Po{v) ~ hiv- v T )F"{v T ) W 

where the primes denote derivatives. The polynomial may be written as 

pm = m + s j c ( , _ f)f * (f) - 2^ i c /(r) La - fjrcfl J r 



" 2* J c ^ F'(f) 



(8) 



where Q(£, v) is a polynomial in v of degree ji — 1, 

and /(i>) has been added to remove the contribution of the pole at 
£ = v. The path C is taken in the counter-clockwise sense and encloses 
£ = v and the zeros of F'{Q but no singularities of /(£)• 
The expression for /(v) obtained from (8) gives 



(10) 



I = J dv i(v) exp [«*•(»)] 

- £ db P (») exp [xF(v)] + f L dv exp [xF(v)] ± £ ^i^p^y 

In order to simplify interchanging the order of integration in the 
double integral, we cut off the tails of L in the usual fashion. The 
error introduced by truncation is exponentially small compared with 
the terms that remain. Deforming C so that it encloses the truncated 
L (in the sense that it encloses the point £ = v for all v's on the trun- 
cated L), interchanging the order of integration, integrating by parts 
with respect to v, neglecting the contributions from the integrated 
portions at the ends of L, and reverting to the original order of inte- 
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gration carries (10) into 

f dv f(v) exp [xF(v)] - f dv P (v) exp [xF(v)] 

Incidentally, if the contributions from the ends (say at a and 6) 
of the truncated L were not neglected, the right side of (11) would 
contain the additional term 



\ \M - PM\ exp [xF(v)] \ 
L xF'(v) _L 



The procedure used to establish (11) can be used to show that, 
for any function /„(£) analytic inside C, we have 

J dv f n (v) exp [xF(v)] 

= f dv P n (v) exp [xF(v)] + - f dv f n+1 (p) exp [xF(v)] (12) 

J l X J L 

where 

UiW 2*iJcF'm-vr' (13) 

Setting / (£) = /(£) and using (12) repeatedly gives 

7=2 x- /" rfy P„(y) exp [atfty)] + R„ , 

^ (14) 



Rk = af 



f>AT 



1 / dvf N+i (v) exp [atfty)] . 



Since Q(£, v) is a polynomial of degree /a — 1 in v, the same is true 
of P n (v) and we write 



P» - YtVnlV , 71 = 0, 1,2, 



i = 



(15) 



1=0 n=0 
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where 

U t (x) = j v l exp [xF(v)] do, I = 0, 1, ■ ■ ■ , n - 1. (16) 

The series (15) is the type of expansion we seek. It would be desir- 
able to have close inequalities for R$, but none are available at the 
present time. 

Another expression for P„(v) may be obtained from (13): 

p aa __If *tQ(t>9) J_ f <fe/-.fe)(-D 

1 n{l) " m io f(r) 2* j c . /?'(«)(« - r) 2 
-artic rffKn L^(r)afJ *"tt) 

In the first line C a must enclose the point z = £ in the 2-plane in addi- 
tion to the zeros of F'(z). Hence initially C s encloses C. When the 
order of integration is interchanged, the only singularity of the inte- 
grand in the f-plane lying outside C is the double pole at £ = z. Ex- 
pand C until it consists of a circle of infinite radius at oo plus a nega- 
tive loop around £ = z. The contribution of the infinite circle vanishes 
because the integrand is a rational function of £ of 0(£~ 3 ) at oo. The 
contribution of the pole at £ = z gives the derivative. 

Notice that the coefficients p„i in (15) arc independent of the path 
L in the v-plane. 

The procedure used to obtain the integral (17) for P n {v) may also 
be used to show that 



'-® - h I df ^~AA) iJUrxr - A ' 



(18) 



When ft = 1 and F(v) = r~, the polynomial P„(v) reduces to 
(- l)7 (2n) (0)/(4 n n!) and j n+l (v) is equal to [/„(») - -vf' n (v) - / n (0)]/(2 V 2 ). 

V. COMPUTATION OF P„[V) , A = 1: BLEISTEIN METHOD 

We shall regard the functions Ui{x) in the series (15) for I as tabu- 
lated or easily computed. For example, when p. = 2 the functions 
Vq (x) and Ui {x) may be expressed in terms of Airy functions. Then 
the most difficult step in applying the series is the calculation of the 
coefficients p nl , I = 0, 1, . . . , /i — 1, of the polynomial P„(v). We 
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desire an expression for p nl in terms of the values of the functions 
g(t), h{t) and their derivatives at the saddle points t = t r ,r = 1, 2, 

. . . , p. 

Let t = t(v) denote the change of variable from t to v, and let the 
saddle point v = v r in the f-plane [F'{v r ) = 0] correspond to t r in the 
£-plane: t r = t{v r ) . We shall use the notation 



K n) = 



An) _ 



:G6VL ■ *• - Kff*L (19) 



'(£) VL • 



r = 1,2, 



When convenient, we shall write h, for /^ 0) = h(t r ) and # r for 0<°> = g(t r ). 
First consider the expression (7) for Po(v). As shown in Appendix 
B, F(v) is a polynomial of degree /* + 1, 

F(iO = £ A iV < (20) 

whose coefficients A, may be expressed as functions of the h r 's. When 
this equation for F(v) is used in (7) , the coefficient of v l in the result- 
ing expression for Pq{v) gives 

Po« = E My Z y//;,, x . « = 0, 1, • • • , m - 1. (21) 

i-l + 2 r = \ r KPt) 

Multiplying the right side of (21) by —1 and changing the limits 
of summarion for ;' from I + 2, p + 1 to 1, I + 1 gives another ex- 
pression for poi- 

Since j(v) = g(t)t ll) , we also need an expression for t { T l) in terms of g(t) 
and h(t). Differentiating F(v) = h(t) twice with respect to v and using 
hp } = leads to 

Kv r ) = g/ r l \ 4? = [F"(v r )/h ( r 2) ]\ (22) 

The sign of the square root is chosen to agree with the constant cmtttcv, 
the form assumed by the change of variable throughout the critical 
region. 

Since the A/s and v r 's may be expressed in terms of the h r 's, equa- 
tions (21) and (22) show that p m depends only on the h r 's, g r 's and 

When n is general, an expression for p nt similar to (21) for p i may 
be obtained by expanding the derivative in the integral (17) for P n {v) 
in partial factions and then using the Cauchy integral theorem. For 
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n = 1 it is found that 

- 2 -,r /(o - /(i>.) f( Pr ) r^l ™ 
'" r L <p. - Vr f + «. - Vr + 2 J (23) 

where the prime on yy denotes that the term for s = r is omitted. The 
primes on f(v) and FCi') denote derivatives with respect to v. 
The expression obtained for P„ (v) is of the form 

PM = t, £ ^ m r- l \v r )/{m - 1)! (24) 

r ■ 1 m - 1 

where a r ( "^ is a polynomial in v. Recurrence relations for the a's may be 
obtained with the help of the partial fraction expansion 



(r - Vry _ ^, (P. - Vr)- [ 1 _ f> (V. - VrY 1 

F'(f) " h f"{ V .) Lr - ». h (r - o ,+1 J ' 



(25) 



The relation ]£' [1/F"(y,)] - - l/F"(v T ) can be used to simplify the 
coefficient of (£ — v,.)-'"- 1 . 
The with derivative of f(v) evaluated at v r , 



i ( 



-«-S(7H(8^i»L (2C) 



contains derivatives of fl(i/). They may be obtained by extending the 
method used to get t$, li . Straightforward differentiation of F(v) = h(t) 
with respect to v leads to 

6t t h T (2?) 

l™ = ijn^m [F w - 6n (8) / (2) / (,)2 - h w f w - dh m t m2 ] r 

where the nth derivative F {n) of F(v) is evaluated at v r and h (n) , 
t in) are evaluated at (,-. 

The values of t u+1) for larger j's may be obtained with the help 
of equation (94) , namely 

F ln \v) = £ h (k) (t)c n . k (28) 

where c n>1 = t (n) , c„„ = i (1,n and the remaining c's are given by the 



1980 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1968 

recurrence relation (96) . Setting k = 1 in (96) gives 

c„ a = n* c -"* cl) + 2 t " 1 )**— , t c - ) (29) 

for w ^ 3, the last summation being omitted when n = 3. The term in (28) 
for k = 1 vanishes when v = v„ t = t r . Substituting for c„, a its value 
given by (29) and solving for i (B-1) leads to the desired result when n ^ 3: 

«-" -spfe" [ F< "' - £ ""' c -'-" <!> § (* » 1 )' , """ , ''" , 1- (30) 

The value of F (n) (v) is zero for ?? > p. + 1. 

VI. COMPUTATION OF P„(v), A = 1: URSELL METHOD 

The Ursell method avoids the evaluation of the derivatives of f(v) 
which appear in equation (24) for P„(v). Instead, it makes use of 
classical saddle point expansions about the individual saddle points 
in the t and v planes. 

Let n different paths of integration, L[, L' 2 , ... , L£ be chosen in (6) 
such that the chief contributions (as x — > °o ) along the paths correspond- 
ing to r, namely L' r in the i-plane and its mate L T in the v-plane, occur at 
the saddle points t = t r and v = v T , respectively. Let the classical 
asymptotic expansions around t T and v r be 

I r = f g(t) exp [xh(t)] dt ~ exp [xh(t r )] S «r»af " H (31) 

[U t (x)] r = f v l exp [xF(v)] dv ~ exp [xF(v r )] £ fti.*""" 1 . (32) 

Using /i(£ r ) = F(v r ), substituting (31) and (32) in the uniform asymp- 
totic expansion (15) with N = «>, and equating coefficients of af*"* gives 

,=0 m =0 (33) 

= TtPrWPn, + E E/5rI«P«-«.« 
1=0 1-0 «•-! 

where the second sum in the last line is omitted when n = 0. The 
expression (17) for P n (v) shows that P n (v) remains the same, irrespective 
of the path of integration L, as long as j(v) and F(v), that is, g(t) and h(t), 
remain the same. Hence, for n = and r = 1, 2, . . . , ju, (33) furnishes n 
simultaneous linear equations which may be solved for p i, I = 
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0,1, ■ • • , ju — 1. Similarly when n = 1, (33) determines the pu's, and so 
on. It turns out (see equation 39) that fi rl0 = w'/3 r0 o . This allows us to 
write the simultaneous equations in the form 

P (v r ) = a r0 //3 r0 n , r = 1, 2, • • • , fi , 34 s 

p.w = f- - t £ (£*)*-.. ■ 

PrOO m = l i=0 \Pr00' 

Expressions for «,-„ and /?,.;,„ may be obtained from the classical 
saddle point asymptotic expansion (103) given in Appendix D. Chang- 
ing n to j in order to agree with the notation of Appendix D gives 

where 7^ n) , (7* B> are the derivatives defined in equations (19), and 
(x) — 1, (a:)„ = xix + 1) ... (a: + n — 1). The 6 m „'s are computed from 
the recurrence relation (100), namely 

^ n-m+l 

Om + l.n + l = ~ i 7 2-1 k a kb m .n-k+l > 

starting with 6 o = 1 and using 

r^U + 2) 

«*= -(/ c + 2) I/,.-'" *-l,2,.... (36) 

The value of arg [— 2//i< 2) ] J is equal to arg {t — t r ) on the portion of L' r 
(deformed into a path of steepest descent through t T ) leaving t r . 
Similarly, 



ft«, - W [_ F < 2) J 2. (2 . _ n)I 2. 6„,„( 2 ), n+ , 

where now the 6„ m 's are computed from (100) with 



(37) 



°< = (*+w- F ' = .(l) FW • (38) 

Setting ;' = in (35) and (37) gives 

" r " = (7r) {^]^°'' ft » = W{^]' y ' > 

where i' 1 ' and /(v r ) are the same as in equation (22). The relation 



(39) 
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P (y r ) = j(v r ) is the starting point for the Lagrange interpolation formula 
(7) in the Bleistein method. 
Setting n = 1 in (34) and ; = 1 in (35) and (37) leads to 

PrOO 1=0 PrOO 

/Sr.o ~ L «" _ I 4 4/> (2) + 9 Ll6/ i (2, + 48/* (2,2 J/ r (4U} 

^"L^Jl 4 " ~ 4F <2) + "Ll6F (2)+ 48F (2,2 J/ r 

where the subscript r on the braces indicates that the enclosed g's, 
h's, v's, F's have the subscript r. 

VII. THE BLEISTEIN METHOD FOR \^ 1 

Here we deal with 

J = f < X ~V(0 exp [xh(t)] dt = J v x ~ l i(v) exp [xF(v)] dv 

Jl- Jl (41) 

M = g(t)(t/v) x -H a) , l (x) = dt/dv. 

The origin is now a singularity, and its vicinity may contribute to J just 
as the vicinities of the saddle points do. Accordingly, we now require that 
the polynomial P (v) be such that P o (0) = /(0) in addition to P (v r ) = 
/(y r ), r — 1, 2, . . . , /*. Assume for the moment that F'(0) ^ 0,that is, 
that the origin is not a saddle point. Starting with Lagrange's interpola- 
tion formula and proceeding as in Section IV gives 

_ /(0)t>F'(iQ f* 1{Vr)vF>(v) 
n{v) ~ vF'(0) + h (0 - » r W"(M 

" /w + z*i i c (v - r)fF'(r) 

%a J c fF'(f) 



(42) 



where C encloses £ = v, £ = 0, £ = ti r , r = 1, 2, . . . , /* but no singu- 
larities of /(£). Here P (v) and Q(£, v) are polynomials of degree 
/a instead of ^ — 1. 






(44) 
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When the origin is a saddle point, P (v) is still given by the expressions 
in (42) which contain integrals. In fact, we have P (0) = /'(0) (with 
primes denoting derivatives) in addition to P o (0) = /(0). 

Much as in Section IV, we obtain 

J tbtP~ l fJp) exp[xF(v)] 

= f dv v x ~ l PJp) exp [xF(v)] + - f do v*- 1 j n+1 (v) exp [xF(v)] (43) 
J i x J i 

where f (v) = f(v) and 

UlW - 2iri J c fF'(f) (r - .) 2 ' 

f - w - a*t Jc r*"(r) ' 

Equations (43) and (44) lead to the desired series for J: 

N - 

•/ = Z af " / ^ " x_1 P»(«) exp [xF(v)] + # A . , 

n-0 J L 

R„ = af*-» f etoi^" 1 / jV+l (y) exp [«F(i»)] dt>. 
When P„ (i>) is written out we get 

Pn{v) = EM', 71 = 0, 1,2, -•• 

J = E PiW E P.i*~" + K* , (46) 

i = n-0 

F,(r) = ( v t+x ~ l exp [.r/<»] dv, J = 0, 1, 2, ■ ■ ■ , M . 
Furthermore, the recurrence relation (44) for f n (v) leads to 

f x ± C X QU, v) 



(45) 



(47) 
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When A = 1 the formulas of this section do not reduce to those of 
Section IV since they contain the additional condition P„(0) = /(0). 
However, (46) gives the same series for J as (14) does for I because 
Vp(x) can now be expressed as a linear combination of V n (x), . . . , 
Vn-i (x) [which become U (x) ,... , £/ A -i (x) ] . 

The only singularities enclosed by C in the integral (47) for P„(v) 
are poles at £ = and at £ = v r , r = 1, 2, . . . , /*. Evaluating the 
integral by Cauchy's theorem gives the coefficients in P n (v) as the 
sum of derivatives of f(v) at v = and at v = v r . The derivatives 
at the saddle points may be obtained by differentiating 

In j(v) = In g{t) + (X - 1) In - + In t w (48) 

V 

with respect to v and using the expressions for #* ] developed in Section V. 
The derivatives f n) (0) may be computed with the help of the series 

t/v = «i" + fj /r + ^ i n + • • • (49) 

where t l n) denotes the nth derivative of t with respect to v at v = 0. The 
C's may be obtained by differentiating F(v) = h(t) repeatedly with 
respect to v and then setting v = 0. If F'(0) =^ 0, 

/l ° (50) 

F 2) - ftTti"' 

where the subscript refers to t = when it is on ft and to v = 
when it is on F. Higher order derivatives may be computed by using 
the results of Appendix C in much the same way as in Section V. 
It may be verified that 

/(0) = g to ^ 

,(i)/ ft \ ._ /<i>M „ ( "/ (1> 4- fo + 1) do ^o 

/ (U) — t g„ to -+- ^ j (1) 

where <7o n) is the nth derivative of g(t) with respect to t evaluated at t = 0. 

VI 11. THE URSELL METHOD FOR \ 5^ 1 

When the origin is not a saddle point, the /i + 1 linear equations 
to be solved for the coefficients p nI , I = 0, 1, . . . , /* in P„(v) turn out 
to be 



an 



'■II 7. ( 1 ) 



/ <2) = 
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PM =f-EE (Hp--' ■ <• = L 2 . • • ■ • " W 

PrOO m-1 1-0 \PrOO/ 

p.o = f=- - 1 s (f4>-... < 53 > 

POUO = 1 J-l-ro-« \PuOO' 

where the summations are omitted when n = 0. The summation condi- 
tion I + m = g in (53) is also subject to ^ 2 ^ m, ^ m ^ » . 

Equations (52) are given by the analysis of Section VI for the 
case A = 1 when g{t) is replaced by t x_1 fif(t), «' by w I + x-1 l /,• by J r , 
and l/|(a;) by V t (x). The «'s and /3's in the rth equation of (52) are 
the coefficients in the classical saddle point expansions about t r and v r : 



-n-J 



(54) 



-r„-i 



,/ r = / t^git) exp [a*(0] ^ ~ ex P [xh(Q] Z <We " 
[V,(a:)] r = f »' +x_1 exp [xF(»)] rf» - exp [xF(» r )] £ fti«a 

J L , m-0 

The a r| in (52) (with j for ?i) is given by equation (35) for a Ti with 
gW-*i replaced by §™~ n) where 



(55) 



*-o W 



The /8 ry in (52) is given by equation (37) for p rU with I replaced by 
I + X — 1 on the right side. 

Equation (53) arises from a consideration of the region around the 
singularity at the origin. As described in connection with equation (106) 
in Appendix D, let L' be a loop enclosing the branch cut running out 
from t = 0, and let L be its mate in the y-plane. Then, as .t- — * «, the 
a's and /3's in (53) are defined by 



Jo = f ' X ~V0 exp [xh(t)] dt ~ £ a 0n x— 

J !.<)■ n-0 



[V,(*)]o = f i 1 **" 1 exp [xF(p)] dv ~ £ ft,,, 

•'lo "'=0 



(56) 






Substituting (56) in the uniform asymptotic expansion for J given 
by (46) and equating coefficients of x~ " _x gives (53) . 

Using the asymptotic series (.106) to determine oo„ and J3 r im leads to 



<*o. 

Pooo 



= '"lWS(f^)iS Mx) "*' (57) 
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where the subscript refers to the origin and b„,„ is computed from 
(100) with a k- given by 

a k = -^ +1) /P+1)!C]. (58) 

The value of # n is the value of dt/dv at v = determined by the change 
of variable from t to v. Similarly, 

PoOO \?0 I m=0 

where, replacing j by n, b„ m is computed from (100) with 

a k = -Fr"/K* + I)' *o"L * ft 1- (60) 

In Appendix F the theory which has just been developed is applied to 
the case of one saddle point (p = 1, X ^ 1). 

So far in this section it has been assumed that the origin is not a saddle 
point. Now let the /*th saddle point coincide with the origin so that F l l) 
and h ( l) vanish. The n + 1 equations determining p nU I = 0, 1, . . . , p 
are now 

«r„ ^ V- (fir 



PM -f-ZZ [T >"-' r - 1, 2, • • • , m - 1 (61) 

PrOO m-1 1=0 \PrOO' 

p..« =f-EE (ffW- (62) 



P„l = o ao.2n+l — /3 0lP n — 2 2^ PolmPn- Q .l 

POlO L O-l m + Z-2«+l -1 



(63) 



where the summations are omitted when n = 0. The values of I and m 
occurring in the inner summations in (62) and (63) must also satisfy 
^ I ^ n and ^ ra ^ » . 

Equations (61) are the same as (52) except that r runs from 1 to 
ix - 1 instead of from 1 to p. The «'s and 0's in (62) and (63) are the 
coefficients in the asymptotic expansions 



J = [ t^git) exp [xh(i)] dt ~ £ a ,af " +X)/2 



[Vi(*)]o = f ^ +x_1 exp [aF(»)] <fo - 2 OI , 

Jrn ».-0 



(64) 



-(m+l+X)/2 



where, as discussed in connection with equation (109), the paths L' , L 
coincide with the paths of steepest descent through t = 0, v = except 
for indentations at those points. 
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Equation (62) is obtained by substituting (64) in the uniform asymp- 
totic series for </<, given by (46) and equating coefficients of aT <2n+M/ \ 
Equating coefficients of x - <8fl+1+X)/a gives equation (63). 

Some results for the case of two saddle points, one of which is at 
the origin, are stated in Appendix H. 

IX. SIMPLE POLE AT THE ORIGIN 

When there is a simple pole at t = and one saddle point in the 
critical region, a case discussed briefly by Bleistein, 3 we have 

J = I T 1 exp [xh(t)] dt. (65) 

In the critical region, U is assumed to coincide with the linear path 
running from a — ioo to a + loo , « > 0. 

Let h(t) be real when t is real and in the critical region. Let the saddle 
point tt. lie on the real axis. As usual, h„ = 0, /if" = 0; and we assume 
/j, g 0, h\ 2 ' > 0. As suggested by example (i) of Appendix 3, we choose 

F(v) = v 1 — 2y,v, v\ = —hi 

where v x is real. We write 

fi =T7 r (-ft,)*, (-/',)' ^ 
I f i I 

in order to make Vi and t have the same sign. 

Equation (46) shows that the uniform asymptotic expansion for 

J has the form 

•/ ~ V Q (x) 2 p B0 .Tf" + VA.v) £p..i"" (66) 

where, with L parallel to, and to the right of, the imaginary u-axis, 

V (x) = [ v~ l exp [xF(p)] dv = ix[l - erf (v^)] 

r,(.r) = / exp [xF(r)\ dv = iU/x)* exp (-xv\) 

erf (z) = -t [ exp (-/ 2 ) dt. 
Putting A = Oin the integral (47) forP„(v) gives 

PM = PnO + PnlV 
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" 2vi Jc 2"f \r - u, ar/ V 



r-* (67) 



/(„) = vl a) /l, i (,) = <B/d» 

where C encloses £ = and £ = Ui but no singularities of /(£). Setting 
v = in (67) gives 

Po(0) = Poo = /(0) = 1, 

P-(0) = Pno = 0, n > 0. 

Here the series (49) for t/v has been used to show that /(0) = 1. 
Therefore the series for J reduces to 

J ~ iw{l - erf [0x3*]} + t"(x/a:) 1 exp (-a»i) £ p»iaf ". (68) 

71 = 

Setting v = v t in (67) gives 



Po, = WW - i]/«. = T" - ,7 



p J_ f fijffi (_J_ ±Y(-X-Y n > 0. 



(69) 



From (22) and F"(iO = 2 it follows that J'" = [2/h[ 2) ] h . The integral 
for p n i may be evaluated in terms of the 2wth derivative of j(v)/v = 
t ll) /t = (d/dv) In t(v). Thus, writing f/(f - t;,) as 1 + y,(f - v,) _1 
and using 

f_ L- AV (r _ Pl )-' = (-i)"i-3 ... (2n - l)(r - v l y 2n '\ n > 
leads to 



(-my 1 

Pnl - y 2„ + l ^ (2rl) . 



7ir £ii _ 

AW « J.-., 



The first of the p„i's required in the series (68) for J is given by 
equation (69) for p i • The remaining ones may be obtained by using 
the Ursell method equation (52) . Since p 00 = 1 and p„ = for n > 0, 
equation (52) gives for n > 



flPlOO L m-l _l 



(70) 
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From equation (55), §(t) = t~ l and its nth derivative at U is §\ n) = 
(-l) n n!<r n-1 . Replacing g by § in (35) leads to 

«., « w(^T +i t (-DV*- 1 £ M*)-, 

L"l J n = m-0 

where arg [— 2/h[ 2) f = tt/2 and the b mn 's are computed from (100) with 
a k - -2/ l , (i+2, /[(fc + 2)! ««], fc - 1, 2, • ■• . (71) 

Equation (54) shows that the /Sum's are the coefficients in the asymp- 
totic expansion of [Vi(.t)]i, that is, in the asymptotic expansion of an 
integral, which has the same integrand as Vi{x), taken along the 
path of steepest descent through the saddle point v = V\. Instead of 
obtaining the /?'s by the general procedure outlined in Section VIII, 
we notice that the "asymptotic" series for V% {x) consists of only one 
term. Consequently fi llm is when m > and the summation in 
equation (70) for p„ x disappears. Moreover, the asymptotic series for 
the error function gives, when v, > 0, 

FoCr) ~ *(»)* exp (-xv 2 ) £ (-l) m (^) m W m_i . 

When v t > and x -» oo , V {x) is given asymptotically by the con- 
tribution from V\. When v 1 < 0, the asymptotic expression for V (x) 
contains the constant term %ri, but the contribution from a path of 
steepest descent through v\ is still given by the same expression as 
for positive V\. Hence, irrespective of the sign of v 1} 

a*. = •Gr) l (-ir(*)^r i "" 1 . 

These results enable us to write equation (70) (with j for n) as 

Pi, = i =^ h IP-'r' ' t (-«.)" t M* + *). - l} (72) 
for ;' > 0. Here, to repeat, 



„ = -L- (- h \* t<» = ixt ^ = — r^i 

Vi k r /iij • ' i uwJ ' i, i *, i L h™ m 

(c) = 1, (c)„ = c(c-f- 1) ••• (c+n - 1) 



(73) 



and the 6„ m 's are computed in succession from equation (100) with 
a k given by (71). 
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The first two p„i's are 

lM" ,1 

P01 = * ~ . (74) 

and the 6„,„'s for p 2 , may be read from the table in Appendix D. 

When | U | is small, the expressions which have been given for p nl are 
essentially small differences between large numbers. If the calculation is 
being performed on a digital computer it may be advisable to use double 
precision. Expanding h(t) about t = I, and then setting t = leads to a 
series for — ft, which may be used to obtain u^ 2 ' -1 as [t[ u /tj] 2 ' +1 times a 
power series in i x . Series of this type can be used to show that 

2 (75) 

Pti = -it[ m $<h + ¥fliOi + S M) + 0(«.) 

where a k is given by (71) . 

X. THE NONCENTRAL x 2 DISTRIBUTION 

Let x be a positive integer and y u y 2 , ... , y s be independent gaussian 
random variables with unit variances and respective mean values 
ft j fti • • • j Z7* • Let ^ De the (noncentral x 2 ) random variable 

= - £ tt • ( 76 ) 

It may be shown that the mean value of z is i = 1 + ?* and that its 
variance is (2 + 4r)/.r where 

•C n-l 

Furthermore, from Rcf. 2 (with a change of variable) the distribution 
function of z is 

Prob [0 ^ z ^ s] = ~ f r 1 exp [a*(0] di (77) 

where d > and 

fe(0 = h\st - ln(l + t) + r(l + 0" 1 - '•]• (78) 
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The integral on the right side of (77) is seen to be equivalent to 
minus Marcum's Q-function (and also to an expression given by 
R. A. Fisher) when (77) is written as 

Prob [0 ^ z <, s] = % [ (z/r) u/i) -> exp [-*(* + r)/2I7 {( ^ l _ 1 [«W)*]A. 

Here I denotes a Bessel function with imaginary argument. 

We are interested in computing the distribution of z when x is large. 
The equation h {1) (t) = gives two saddle points. However, as pointed 
out in Ref. 2, when x is large only the one at 

1 + (1 + 4rs) h 
h ~ 1 + 2s 

need be considered. The value of t t is real and > — 1. When s = 1 
+ r, t x is zero; and when s increases through 1 + r, t t decreases 
through 0. 

From equation (68) the desired asymptotic expansion is 

Prob [0 ^ z ^ s] ~ |{1 - erf [v^]} 

+ !(«)-* exp {~xv\) £ p nl x~ n (79) 

n = 

where p nl is given by (72) . The quantities entering p nl are 

1 4 i i 

fc w (0 = K-l)'(n - 1)! [nr(l + i)""" 1 + (1 + 0~"], n ^ 2 

M" = ki + «.r 3 + 2- J (l + g- 2 



2(- 

a, = 



l) t+1 ft* + 2)r + 1 + ft ] , > , 

(1 + 0*1- 2r+l + «i J' " ' 



(ft + 2)(1 

The values of p 01 and p u may be obtained from (74) by substitution 
of the parameter values (80). 

When s — z = s — 1 — r is small, the central limit theorem in the 
theory of probability states that 

[" (g - z)x i 
_(4 + 8r)*J, 

This agrees with the approximation given by the error function term in 
(79) when it is noted that t t « — (s — 1 — r)/(l + 2r) and — /i x « 
<?A{ 2) /2 if s - 1 - r is smaU. 



Prob [0 ^ 2 ^ s] ~ | i 1 + err 
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The ordinary x 2 distribution is obtained by setting r — in the 
noncentral x 2 distribution. In this case we have 

«i - ^ , *r = 2/ S) 

o 

y ' = |T=7| [ 2 J ' (81) 

2(-s) fc . . , 

flfc = -fcT2' fc = L 

Equations (74) show that the first two coefficients p i, Pn in the 
asymptotic series (79) are now 



2 1 

1 U 2 



Pu = 



:-^]-i}- 



(82) 



2 IA1 - s; 
When s is close to its average value 1, equation (75) gives 

pox - 2/3 + 0(0 
p„ = 1/135 + 0(0- 
Setting a; = 2c and r = gives 

1 /* e * 
Prob [0 ^ z :g s] = -=tt / w c * exp (—it) tfi* 
l {c) Jq 

_ Y^ (cs) n exp (— cs) 
" Ai n! 



(83) 



(84) 



where c is assumed to be an integer (x even) in the last equation. 
These relations may be combined with the foregoing formulas to 
obtain asymptotic results for the incomplete gamma function and 
the Poisson distribution. 

There is reason to believe that the asymptotic expansion (79) for 
Prob [0 ^ z ^ s] may hold over the entire range fS s ^ °o . For exam- 
ple, consider the ordinary (r = 0) x distribution. In this case the first 
two terms in (79) give 

Prob [0 ^ z ^ s] ~ |{1 - erf fax*}} 

exp (-xv 2 ) (85) 



+ *[«]"* [j 2 l 



— S V 



1- 
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where x = 2c and v Y is given by (SI). Let c be held fixed at some large 
value and consider further the three cases s — » 0, s — > °o, and s — > 1. In 
all three cases it may be verified that the leading terms in Prob 
[0 ^ z ^ s] given by (85) agree with those obtained from the exact 
equations (84) and the asymptotic properties of the incomplete gamma 
function. The expressions obtained from (84) are 

Prob [0 ^ z S s] = (27rc)" 5 s r [exp c] [1 + 0(cs) + 0(0], s -> 
Prob [0g«|8] = l- (27Tc)-V" , [exp(c-cs)][l + 0(s-')+0(c~ 1 )], 

s -> oo (86) 

Prob [0 ^ e ^ 1] = § + |(2wc) H + 0(0, s - 1. 

APPENDIX A 

77ze Behavior of h(t) in the Critical Region 

In this appendix we show that, in the critical region, h{t) behaves 
much like a polynomial of degree /x + 1, and we examine the change 
of variable from t to v. 

First write 

/,(/) - £ i^j- + fl„ + , (87) 

where /t,'/' is the value of (d/cltyh(t) at « = 0, /i< M+,) is not 0, and R, +l 

is 0(0- 

One of the distinguishing features of a polynomial in t of degree 
n + 1 is that when I is much larger than r, where \t\ = r is the smallest 
circle which encloses the zeros, the dominant term in the polynomial is 
the one containing t" + \ The function h(t) has a corresponding property. 
Suppose that the saddle points tj,t 2 , ■ ■ ■ , t„ all he within a distance e of 
the origin, and for the moment suppose that they may be moved towards 
the origin so that e may be made as small as we desire. Also suppose that 
fcj* ,+1> = A + O(e) where i^0. We shall show that by making e small 
enough we may find a range p < | t \ < 77 throughout which 

h{t) » «" + % +1 /(m + 1)!. (88) 

Here 77 is such that when 1 1 | and e are less than 77, the remainder R „ +1 in 
(87) is negligible in comparison with the (p. + 1) term T „ +l = t" +1 ho li * 1) / 
(n + 1) !. Once 77 is fixed, p may be chosen to be arbitrarily small, subject 
only to p < 77. 

In order to show that (88) holds when p < | t \ < 77, notice that by 
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repeated differentiation of the representation 

h'(t) = (t - t t )(t - t 2 ) ... (t- t,)G(t) 
it may be shown that 6(0) = K* +1) /ii\ + 0(e) and that h« 
0[e* +w ^ M+1) ] for j = 1, 2, • • • , ii. Hence when \t\ > e, 

£ M!. = [t" e C +1) ] - 0(7 r„ +1 



IS 



1-1 



Choosing p to be arbitrarily small, subject only to < p < v, and then 
choosing e so that e/p « 1 establishes (88). 

This property of h(t) suggests that some insight into the change of 
variable from t to v, specified by F(v) = h(t), may be obtained by 
considering h(t) to be a polynomial <p(l) of degree p + 1. Then a natural 
choice of F(v) is F(v) = <p(cv + 6) where c and 6 are constants. For 
simplicity we take c = 1 and b = so that w r = t T , r = 1, 2, . . . , n 
where v, is the rth saddle point on the u-plane. The equation F(v) = h{t) 
goes into <p{v) = <p{t) which we write as 



v(v) - <p(t) = £ A& - n 

;" = 1 

= (v - t)[A x + A 2 (v + t) + A,(v 2 + i>J + + • • •] 

= 0. 

The branch used in the change of variable is 

v = t (89) 

for which dv/dt = 1 everywhere. The remaining /a branches, which 
are ignored in the change of variable, may be obtained by solving 

.4, + A 2 (v + t) + . . . + A, +1 (v" + v»-H + ... + n = (90) 
for v as a function of t. Writing (90) as 

G( V , t) = *> - f > = 

and expanding G(u, i) about »■{,,<■(, shows that, near t = £ r , one 
of the remaining branches behaves like v = t r — (t — t,). On this branch 
v = t r and dv/dt = — 1 at t — t r . Again, let v m d„ s ■■ 1, 2, . . . , n — 1, 
be one of the /x — 1 roots of (?(?;, t r ) = which is not equal to t r . Expand- 
ing G(y, t) about v = i)„t = t T shows that near £ = f r the correspond- 
ing branch behaves like 
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2<p (y.) 

and that dv/dt = at t = t r . 

The examination of special cases suggests that there is a one-to-one 
correspondence between the /t branches of (90) and the /x saddle points 
in the sense that dv/dt for a particular branch is equal to —1 at its 
corresponding saddle point and is zero at the other saddle points. Thus 
it appears that the branch v = t [or its analogue for general h(t)] 
is the only one suitable for the change of variable throughout the 
entire critical region. 

APPENDIX B 

The Choice of F{v) 

The polynomial F{v) used in changing the variable of integration 
will be written as 

F{v) = 2 A iV \ (91) 

,-0 

The positions t\ , t% , . . . , tn of the saddle points and the associated 
values h r = h(t r ) are supposed known. We require expressions for 
the A/s which are either pure numbers or depend only on the h r 's. 
Although one or more of the zeros v lt v 2 , . . . , v„. of F'{v) = dF{v)/dv 
may appear in our final expression for F(v), they will always be ex- 
pressed in terms of the Vs. 

Since F(v) = h (t) , we have the 2//. equations 



F(P r ) = h T , 

F'(v T ) =0, r = 1, 2, 



(92) 



relating the 2/x + 2 unknowns v u v 2 , ... , v „ A Q , A u ... , A M+1 . Con- 
sequently we have at least two arbitrary choices (A ,,+ ! = is forbidden). 
For the case X =^ 1 we shall always require the change of variable to be 
such that v is when t = and thus we take A = 0. In some instances 
the form of h(t) aids in the choice of the A ,-'s. For example, when h(t) is 
an even function of t, we can take F(v) to be an even function of v. 

In choosing F(v) it is helpful to notice that in the critical region the 
change of variable takes the form t tt cv, c being a constant. Conse- 
quently, from (S7), 



"+1 /„ A »'!,<'> 



F(v) » E *?&**_ (93) 
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and it follows that the /x zeros, v r , of F'(v) have nearly the same configura- 
tion in the w-plane (except for a possible rotation and magnification 
given by arg c and \c\ , respectively) as do the zeros, t„ of h'(t). For the 
case X = 1 there may also be a small displacement so that t « cv can be 
written more accurately as t ^ cv + b, or dt/dv ^ c. Furthermore, from 
(88), we can take c to be one of the roots of 

A, +1 =c h+1 ^ +1 7(m + 1)!. 

The following examples illustrate possible choices of F(v). 
(t) n = 1, X 5* 1, h 9^ 0. Initially there are 4 unknowns, v if A , A lt A a 
related by 2 equations, and F(v) given by 

F(v) = A 2 v 2 + A iV + A . 

We take A = (because X ^ 1) and arbitrarily choose A 2 = 1 (for 
convenience). This carries the two equations into 

v\ + A 1 v 1 = h x 

2», + A, = 0. 
Consequently 

ify) = v 2 - 2t>iti v* = -/h. 

This case has been considered by Bleistein. 3 

(ii) n = 2, A = 1. Initially there are six unknowns Vi , Vo , A ,A X ,A 2 ,A 3 

related by four equations, and F(v) given by 

F(v) = A 3 v 3 + A 2 v 2 + A,v + A . 
We take A 2 = in order to simplify F'{v) . The four equations become 

F'(v r ) = SA 3 v 2 T + A, = 

3F(» r ) = -A^ + 3A,y r + 3A = Sh r , r = 1,2. 
It follows that v 2 = -«i, A = (hi + fo)/2 f A x v x = 3(/ii - /i 2 )/4. 
For the remaining choice we take A 3 to be equal to — Ai/3 and 
obtain, 

F{v) = \Qh - 1h)(v 3 - 3v) + h(h 2 + /'.)• 
Another choice for A 3 is ^ which gives 

TO = & ~ v\v + i(A» + *,) 

y? = l(/i 2 - fcOi y 2 = — »i • 

This case has been considered by Chester, Friedman and Ursell. 6 
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(Hi) n = 2, X ^ 1, U t± 0, t 2 = 0. Here the unknowns and F(v) are the 
same as in example (ii), but because of h 2 = h(0) = it turns out we 
have three arbitrary choices. Since X ^ 1 we take A = 0. Then both 
F(v 2 ) = h 2 a.ndF'(v 2 ) = are satisfied by choosing v 2 = and Ax = 0. This 
leaves F(vi) = ^ and F'(ih) = to be satisfied by the remaining three 
unknowns v u A 2 , A 3 . Our third choice is A 3 = 2. It leads to 

F(v) = 2v 3 - 3v/, v\ = -h, . 

(iv) y. = 2, X ^ 1, U, t 2 , t^ 0. This case illustrates the complication 
encountered for the general case when p ^ 2. The value of A must be 
and we choose A 3 = 1. Then 

F(v) = v 3 + A 2 v 2 + A Y v 

F'(v) = 3v 2 + 2A 2 v + A, . 

The last equation shows that A» = - 3(Wi + t> 2 )/2, A x = 3f it' 2 . Sub- 
stituting in F (v r ) = h r gives 

v 2 (- Vl + 3» 2 )/2 = ft, 

vl(-v a + 3v,)/2 = /i 2 . 
Setting a = /i 2 //ii and p = v 2 /vi leads to 

p 8 - 3p 2 + 3pa - a = 
which has the three roots 

p„ = 1 + (1 - a) J [(l - a 5 ) J a; n + (1 + aVco*] 

where w! = 1, a 2 = * 4/3 , w 3 = i~ 4/3 and the star denotes "conjugate 
complex." When U and t 2 tend to zero, one of the p„'s tends to t 2 /t lt and 
this is the value of p to be used. The value of v\ is equal to 2hi/(Zp — 1), 
and we have 

A 2 = -3»,(1 + p)/2, A x = 3pvl 

(v) p. = 3. The general case of three saddle points may be handled 
by a procedure similar to that used in example (iv) . We do not dis- 
cuss this case beyond mentioning that when we set v 2 = pfi , v 3 = aVi 
the variable u = (p— l)/(o— 1) must satisfy the equation 

u 4 - 2u + 2au - a = 0, a = (h t - h 2 )/(h } - h 3 ). 
(vi) p = 3, fa = 0, M0 even. Since M0 is even we start with 

F(v) = v 4 + A 2 v 2 + A n 
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and find that 

F(v) = v 4 - 2v]v 2 , v\ = -h^ 
This case has been treated by Felsen. 7 

appendix c 

Derivatives of Composite Functions 

A result used in Section V to compute the nth derivative, tl"\ of t(v) 
at v = v T is stated in this appendix. Let the argument u in h(u) be a 
function t(v) of v. Then 

(£fh[t(p)] = £ Wi , n fc 1 (94) 

where h w stands for (d/du) k h(u) evaluated at w = t{v), and the 
coefficients c n< k are computed from the recurrence relations 

c n .i = t in \ (95) 

, = T,( n )t ln+1 - m) c m , k , l^k^n (96) 



n 
! ; "" 7=h \m. 



in which t ln) denotes (d/dv) n t(v) and I ) the binomial coefficient. 
Equation (96) may be proved by induction. Differentiating (94) gives 

W*+i = t tl} c..k + |c n . A . +1 , 1 £ * £ n - 1. (97) 

We assume that (96) holds when n is replaced by n-1 and use it to 
express c„, fc+ i as a sum. Then one of the terms in the summand for 
d c Hi k+i/dv contains d c m ,je/dv. From (97), assuming k > 1, 



— r - r — t w r 

Equation (96), with (n— 1, k— 1) for (n, fc), lets us sum the terms 
containing £ (1) c m , fc _i with respect to m. Equation (96) follows upon 
combining binomial coefficients and using c ki k = £ (1, c»-i, *_i. 

The recurrence relations may also be obtained by writing the right 
side of (94) as a Bell polynomial and using the recurrence relation 
for these polynomials. 8 Expressions for the c„, k s (up to n = 8) as 
polynomials in the t {n), s may be obtained from Riordan's table of Bell 
polynomials given on page 49 of Ref. 8. 
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APPENDIX D 



Formulas for Classical Saddle Point Asymptotic Expansions 

A result useful in obtaining asymptotic expansions of integrals is 



f t' _i (?(t) exp [xH(t)] dr 
Jo 



AlMV'" ; 



Z <?,-„ Z b mn T 



m 



P + J 



where x -» oo , Re p > 0, i» is a positive integer, and 

G( T ) - £ r n C n , ff (r) = Z T n tf n . 



(98) 



(99) 



The b mn 's depend only on H (t) and are computed in succession, start- 
ing with 6 o = 1 and b 0n = for n > 1, from 



^ n-m+l 

&m+i.n+i = — TT Z kakb m .„- k +i ■ 

71 -t" 1 A = i 



(100) 



Here a fc = — H v+k /H v , k = 1, 2, . . . . 

Special values of 6„,„ are given in Table I. 

The asymptotic expansion (98) is based upon the gamma function 
integral 



and the expansion 



/ u'~ l exp (— u") du = - r(- 
r co -i 

exp y z a»r = Z f Z &m„?/ m . 

L 1 _ n=0 m=0 



(101) 



(102) 



The recurrence relation (100) may be obtained by differentiating 
(102) with respect to £, replacing the exponential by its series, and 
then equating coefficients of $"y m+1 . 

Table I — Special Values of b mn 



» 


ban 


bin 


b?n 


(l3n 


6ln 


6„„ 





1 












1 





ax 










2 





a* 


Oi*/2 








3 





03 


aia-> 


fli 3 /6 






4 





04 


a,a 3 + 2~ 1 a 2 2 


2- 1 a 1 2 a 2 


OiV24 




71 > 1 





a„ 








ai"/n! 
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For an integral, say J r in equation (31), having a simple saddle point 
at t = t T and a path of integration L' r which runs up to, through, and then 
down from t r , we can use (98) with v = 2, p — 1, t = t — t r , G(t) = 
g(t r + t), and H(t) = h(t r + t) — h(t r ). The contribution of the saddle 
point is obtained by deleting the terms in (98) for which j is odd, 
doubling the terms for j even, and taking arg (— l/xH 2 ) = arg 
[— 2/xh { r 2) ] i to be equal to arg (t — t r ) on the part of L' r just leaving t r . 
The result is 



f g(t) exp [xh(t)] dt 



U 



exp [xh r ] zl 



9 



j + 4 2) „<2/-n) 



2 7# r, Z b mn T(m + ; + £ (103) 



^oLC'. »^(2i-n)!ir 

where the derivatives of g(i), A(t) are denned in equation (19) and 
6„,„ is computed with 

<" = - (t + 2)m» - (104) 

The gamma function T(m + j + ^) may be written as Vtt (|)m+,-- 

For the integral J given by equation (56) most of the contribution 
comes from the region near the branch point at t = 0. When t = is not 
a saddle point, there is only one path of steepest descent {for exp [z/i(£)]} 
leaving t = 0. This path may be taken to be the branch cut in the f-plane 
and the path of integration L' n for J may be taken to be a positive loop 
enclosing the cut. Then the asymptotic series for J may be obtained 
from (98) by setting v = 1, p = X, r = t, G(t) = g(t), H(j) = h(t) and 
using in place of (101) the integral 



r 



exp (-») du = [1 - exp (-i27rz)]r(z) = ^g _ , 



(105) 

Here arg w is on the part of the path of integration leaving t = 0. 
The positive real it-axis in (105) is a branch cut. The path of integra- 
tion starts at u = + oo on the top side of the cut, comes in along the 
cut, encircles u = in the positive direction, then runs out to u = + oo 
along the bottom side of the cut. 
It is found that 
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Jo = [ < x_, ff(0 exp [xh(t)] dl 

J/ '°' (106) 

1 ~1 X+ ' » n (i ~ n) 



■ J2 b mn [l - exp (-27rf\)]r(m + X + j) 

m = 

where arg [— l/xh ( x } ] is equal to arg t on the part of L' leaving t = 0, and 
b„„ is computed with 

a * = -(Jfc + l)!*i l >- (107) 

The last relation in (105) may be used to handle the case in which A 
is or a negative integer. 

When I = is a simple saddle point as well as a branch point, the 
path L' can be taken to coincide with the path of steepest descent 
through t = except for an indentation at t = 0. The indentation is 
chosen so that a man walking in the positive direction along L' would 
have the point t = on his left. We put v = 2, p = X, G(r) = g(t), 
H(t) = h(t) in (98) and use in place of (101) the integral 

f .-i / 2n , in / • \it( z \ iw CX P (~ *W2) 
/ u* exp (—vT) du = |[1 - exp {-tirz)]T\^J = - r — 



^'2 

(108) 

Here K runs from u = — oo to u = + oo with a downward indentation 
at u = 0, and arg w is on the part of K leaving u = 0. Instead of 
(106) we now have 

» x /_q\(X+/)/2 ," /7 ( ' _n > 

• £ & mn |{l - exp [-ix(X + i)]|r(m + ±±l) (109) 

7/1=0 ^ - 1 ' 

where arg [— 2/xh l 2) ] i is equal to arg t on the part of L' leaving t = 
and b m „ is computed with 

2h (k+2) 
a * = "(A + 2)! AT' (U0) 
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APPENDIX E 

Two Saddle Points 

In order to illustrate some of the results of Sections IV, V, and VI, 
we consider the case of two saddle points, p = 2. This case has been 
discussed by Chester, Friedman, and Ursell. 4 - 6 From (15) the desired 
expansion is of the form 



/, 



0(0 exp [xh(t)] dt ~ U {x) £ p n0 x n + Ufa) £ Pm* 



(111) 



PJp) = P„o + PnlV 

where, from example (u) of Appendix B and equation (16), 
U t (x) = J v l exp [xF(v)] dv, I = 0, 1 

F(p) = i v 3 - v\v + A , (112) 

A = \{h 2 + h x ), v\ = l(h 2 - fcOi v 2 = - y i • 
Arg Vi is determined by the correspondence of Vi with t x which comes 
with the change of variable from t to v. 

Let U and the change of variable from t to v be such that L runs 
in from v = oo exp (—171-/3) to the critical region near v = and 
then out to oo exp(iir/3) (it may be necessary to reverse the direc- 
tion of L') . Then 

U (x) = 2irix~ i Ai(x i v]) exp (xA ) , ,-_„» 

C/ifr) = -2irtaf 'iit'fcM) exp (xil ) 
where .Ai(z) is the Airy function and Ai'(z) its derivative with re- 
spect to z. 

From F'(v) = v 2 — v\ and equation (9) it follows that Q(f, v) is f + v. 
Consequently equation (17) for P n (v) gives 

*»-a/.*4Fhs]7^ (114) 

Here /(v) = 0(i)< <0 in which f (1) = (tt/cfe is obtained by differentiating 
F(v) = h(t) with respect to v. The path of integration C encloses f = ±Vi 
but no singularities of /(f) . 

From v 2 = — Vi and 

??oo + v r p . = P*(v r ) = /(»r), r = 1, 2 
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we get 



POO = 9 I VOX 7^ (Hi>; 



/(»>) + /(-Pi) _ /(»,) - /(-Pi) 

2l>, 

which may also be obtained from equation (21) for p i- 

Putting v = f>i, n = 1 in (114) and expanding the integrand in 
partial fractions leads to 

PM = jjjj [/(-.,) - /(<<,) + 2v 1 f\v 1 ) - 2v\f\ Vl )]. (116) 

Expressions for p™ and p n which agree with equation (23) for p lt 
may be obtained from (116) by changing the sign of Vi to get Pi( — i>i) 
and using n ■ = 1 in 

PM + Pn(-V t ) 

VnO = 



2 

P.(W - 

Vnl = 



(117) 

PM ~ Pn(-Vl) 
2V, 

If we were to continue with the Bleistein method we would have to 
evaluate / (n) (v r ) by using equations (26) through (30) . Instead we 
turn to the problem of obtaining P„{v r ) by the Ursell method. For /x. = 
2, equations (34) become 

PM = a rO /0 rOO , r = 1,2 



(118) 



PrOO m-1 PrOO 



where a,, and /? r j ; - are given by equations (35) and (37). 

Since g(t) and h(t) in the original integral (111) are quite general, we 
use equation (35) for a rt as it stands. However, equation (37) for (J Tli 
simplifies considerably. This is to be expected since it gives essentially 
the coefficients in the asymptotic expansions of Ai(z) and Ai'(z). From 
Fl 2) = 2v r , F ( r 3) = 2, and Ff" = for n > 3 we have a x = -l/(3» r ), 
and a k = for fc > 1. It turns out that 6 mn is for m ?* n and b nn = a\/n\. 
When I is set equal to in (37), all terms vanish except the one for 
n = 2j, m = 2j. When I = 1, all terms vanish except those for n = 2j, 
in = 2) and n = 2j — 1, m = 2j — 1. It is found that 

«J 9<(2fl! (n9) 



00 ~ \ Vr ) ' &'«* 





2004 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1968 

Changing P n (v r ) into P k (v r ) in order to avoid confusion with the 
n used in Appendix D carries (118) into 

Here r = 1, 2; v 2 = — V\\ and 6 mn is computed from (100) with a k given 
by equation (36). The value of t™ is given by 

where arg t[ l) is calculated either from (i) the form t « cy assumed by 
the change of variable in the critical region or from (ii) arg [ — 2/h l r 2) ] h 
and arg {-l/v T ) h being equal to arg (t — t r ) and arg (v - v r ), respec- 
tively, on the portions of the paths of steepest descent leaving t r and 
v r . The last summation in (120) is omitted when k = 0. The expression 
for Pi(v r ) may be written with the help of equations (40). 

Equation (120) was checked by using it to obtain the first few terms 
in the known 10 uniform asymptotic expansion for the Bessel function 
H™(xz) with < z < 1. Here h(t) = z sinh t - t, the saddle points 
are at ±«i (U > 0) on the real axis, and the path of integration runs 
from t = — oo to t = oo -f iir. If the direction of the path of integration 
is reversed [so that (111) gives -H l x l> (xz)], the paths V and L can be 
brought into correspondence by a rotation of 120°. In the approximate 
form t « cv of the change of variable, arg c = 2tt/3; and v x corresponding 
to <! is f>i = | -3/i,/2 | J exp (-«'2tt/3). Furthermore, /(») turns out to 
be an even function, (114) gives P n (—v) = (-l) n P n (v), and p 2ntl , 
p 2n+ i, are zero for n = 0, 1, 2, • • • . 

When U and t 2 approach each other, h[ 2) , h l 2 2) , v t and v 2 tend to zero. 
In this case the asymptotic behavior of the integral (111) may be 
determined with the help of the equation obtained by setting v = 3 
in equation (98). However, if one is interested in the behavior of the 
coefficients p„, , the following relations are useful. Putting v, = in the 
integral (114) for P„(v) shows that in the limit 

Poo = 1(0), Vox = /'(0), (122) 

p,o = -/ (3, (0)/3! , p„ = -2/ (4, (0)/4! 
and so on. The derivatives t [n) appearing in the derivatives of f(v) = 
g(i)i (1) are now obtained by differentiating 3 _1 v a = h(t) repeatedly 
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with respect to v. The leading coefficients are found to be 
Poo = 0i '1 h — Y^/hi J 

„ Cl).(l)2 I „<0) v (2) .(2) L(4),[l)S/in 

Poi = 0i ii + ffi 'i 'i = — hi ti /12. 



(123) 



APPENDIX F 

Saddle Point Near Branch Point 

Here we apply the theory of Sections VII and VIII to a case discussed 
by Bleistein, 3 namely, X ^ 1 and n = 1. The paths V , L and the func- 
tions h(t), F(v) are assumed to be the same as in Section IX where the 
singularity at the origin was a simple pole instead of a branch point. 
In the critical region L' is parallel to, and to the right of, the imaginary 
£-axis. Only the case in which t Y and v x are real and of the same sign 
will be considered. When U and i\ are positive the cut associated with 
the branch point is assumed to start out from the origin along the posi- 
tive real axis, and then quickly bend downward to run out to —i<x>. 
When U and v x are negative, the cut starts out along the negative axis 
and then bends downward to — i"oo. 

Equation (46) and F(v) = v" — 2v r v lead to 

J ~ V (x) E p n0 x-° + V l (x) £ p ni x~ n (124) 

n=U n=0 

where, with c > and the path of integration lying to the right of 
the cut, 

2 ^ (-2v,x h y 



nC+ ICC 

V (x) = y x_1 exp [xF(v)] dv = u 

•'c-ioo 



,:,U-^± n 



Replacing A by A + 1 gives V\(x). V (x) and Vxix) are parabolic 
cylinder functions (Bleistein, 3 and pair No. 740.2 in Campbell and 
Foster Table") . 

PM = Pn0 + Pni v (125) 

dt /(f)f X_1 



= -• f 

2iri Jc 



2" 



1 



n (r + v - v,)r x 



where C encloses £ = and £ = V\. Setting n = and using f[v) = 
g(t) (t/v) x ~n^ leads to 

Po(0) - Poo = /(0) = 0o O) /o m , C = ~2vM l) 

pm = m> t™ = [2Ai 2> ] i . 
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Setting n = 1 gives 

Pi(o) = ^ [/(o) - m + Vl rm 

PM = ^ [/(0) - /(*,) + v 1 f l \v 1 )) - ^. 

The values of p n0 and p ni may be obtained by following the steps 
outlined in the first part of Section VIII. Equation (52) and (53) 
become 



Pnfa) = /- Un - £ 0S 1O ,«P n - m .O + PllmPn-m.l) 

PlOO L m = l -J 

PnO = o «0n — 2 (/3 0«Pn-,.0 + A, . , ,,. ,£„-«, , i) . 

Pooo L 9=1 - 1 



(126) 



The path L x in (54) is parallel to the imaginary axis and passes 
through the saddle point at v = v x . The path L in (56) runs up along 
the right side of the cut, encircles the origin in a positive direction, and 
then runs down to v = — ioo along the left side of the cut. 

To get fiuj we replace I on the right side of (37) by I + a — 1 and 
notice that a k given by (38) is for the values of k used in computing 
b mn . Consequently, the only nonvanishing b„,„ is b 00 = 1 and (37) 
gives 

PlOO J* 

The expression for «,„ obtained by replacing g\ 2 '~ n) in (35) by §[ 2 '~ n \ 
where §(t) = t x ~ l g(t), contains [— 2//i, (2) ] i which may be written as U\ i} . 
In equation (57) for a 0j //3ooo we replace b mn by S mB to indicate that 
5 m „ is computed with a k given by (58) instead of the a k used in computing 
a ln //3ioo • A still different a k , given by (60), is used to compute /3 i,//3ooo • 
From equation (60) all of the a k 's used to compute p" j,//3ooo are zero 
except cii = l/iZVi). Therefore b mn is zero unless m = n, and it follows 
from (59) that 

/W0ooo = (x)i + 2f(2» l )~ , - a 7;i. 

When all of these results are used in the expression (126) for p n0 we 
get, with ;' for n, 

Vio = (2 yi r^ <x+n Eff^j t L..W-, 

n = \J ~ n )' m-0 

- £ £t& [(X + 2s - l)p/-..o + 2^,-..,] (127) 
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where /i(0) = h = 0, h, < 0, h™ > 0, i\ = U \ U | _1 (-/«i )*, «5" = 
—2Vi/ho\ and &«n is computed from (100) with a k replaced by d k where 

When j = the summation with respect to s is omitted. Similarly, 
P»(v\) gives 

Vn = (-D "■ U 2- (2i _ „)! £ 6 -»<*>-+' y, 

-£ tJ SteM = ¥ l! —-»*~] (128) 

where if 1 ' = [2/li ( r'\% §(t) = « x "V(0i and 6 mn is computed from (100) 
with 

a k = -2hl k * 2) /[{k + 2)\h[ 2) ]. 

It is interesting to notice that when \ = 0, (127) and (128) give 
values of p ;0 and p n which agree with those obtained in Section IX. 

When the saddle point approaches the branch point, t { , h x , h ( l) , and 
v x tend to zero. In this case the integral J may be evaluated with the 
help of equation (109). The behavior of the coefficient p nl may be 
studied by putting v 1 = in the integral (125) for P n (v). It is found that 

Pu„ = /(0), Poi = / U) (0), (129) 

Pu> = -M (2> (0)/2, Pll = -(X + l)/ (3> (0)/12 

and so on. The derivatives t" n appearing in the derivatives of f(v) = 
g{t) (t/v) x ~H il] are now obtained by differentiating v- = h{t) re- 
peatedly with respect to v. 
For n = 1 and 2, 

t!° = P/A? 1 ]*, t™ = -JW/6. (130) 

Substituting these values in equations (511 for /(0) and / (1) (0) and 
using (129) gives the limiting expressions for p 0Q and p 0l . 



APPENDIX G 



Poisson-Charlier Polynomial 

In this appendix equations (1271 and (128) for p j0 and p n arc used 
to obtain an asymptotic series for the Poisson-Charlier polynomial 
c,Ay, a) when y is 0(1) and both n and o are large and positive. 
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Multiplying the generating function 

(1 - u)" exp (au) = £ ~~ C m (y, a) (131) 

by vr"- 1 and integrating u around a small circle enclosing u = gives 
a contour integral for c n (y, a). Instead of c n {y, a), we find it more 
convenient to deal with the polynomial (in y) 

d.(y,a) = a " eX P{- a) c(.V,a). (132) 

lb • 

Setting x = n + 1 and making the change of variable t = 1 — u in 
the integral for c n (y, a) leads to 



1 /*+" £ exp(-aQ 

^' a) = 2wJ WM (i-o* *■ 

= -i, f P exp [.r/i(0] <ft 

2x1 Jl-ioo 



X > V + 1 



2iri ' 



(133) 



fc(«) = -ri - In (1 - 
where r = a/.T. The ratio r is positive. 

We wish to use equations (127) and (128) to compute p„ , p»i in 
the expansion obtained by dividing (124) by 2tti, namely 

dn(y, a) = J~ V (x) £ p BO af " + fife) Z p„ia;"". (134) 

n = n = 

Here f (x) = V (x)/2iri with X = 7/ + 1. The function fi{x) is obtained 
from P («) by increasing y by 1. 
We have 



PoOr) = A (' '" *>" exp [.r(y 2 - 2v,v)] dv, c> 

Z7Tt Jc-iao 






(?[y, ^ 



(135) 



where arg V is at V = c and G(i/, z) is the parabolic cylinder func- 
tion 

C(2/> z) = —. \ u" exp (w 2 — 2zu) du 



= 2-Z 



(-*)" 



«!lH^P 
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The function G is related to the function U discussed and tabulated 
in chapter 19 of Ref. 10 by the equation 

9 -i//2 

G(y, z) = ^t-tx exp (-K) U[-y - |, *2*]. 

The saddle point t = f., is obtained by setting the derivative h ll) (t) = 
_ r _|_ (i _ ty 1 to zero; and the saddle point v = w, is given by the 
relation v] = — /h together with the condition that u, and t x be of the 
same sign: 

t x = \--, vy = rn<r- 1-lnr)*. (137) 

The values of the derivatives t (1) = dt/dv at v = and v = v x are 



,(i) 2»i 



r - 1 

*!" = P/A?']* = ^ 

For k > 1 the fcth derivative h ik) (t) is (fc - 1)!(1 - t)~ k and the 
coefficients used to compute 6 m , n , 6 m . n from (100) are 

d k = -k\/[(k + 1)! (1 - r)] = 1/P + l)(r - 1)] 

and 

a, = -2(Jb + l)!^ +i /[(* + 2)!r 2 ] = -2i^/(* + 2 )> 

respectively. 

Comparison of the integral (133) for d.(y, a) with the integral (1) 
for J shows that g(t) = 1 and X = y + 1. Consequently £(<) = £ x l flf(0 
becomes §(t) = t". For k > the derivatives are g (k) = and § (k) = 

y(y - 1) • • • <y - * + «'"'• 

Setting ;' = in (127) and (128), and using the results just ob- 
tained gives 

( 2 ''- Y +1 fcV@2 Po ° mw 

for the leading coefficients in the asymptotic expansion for d n {y, a). 
Here t x and v t are given by (137). The next two coefficients, obtained 
by setting ; = 1, reduce to 

p.u - (y + 1)0/ + 2)Poo|_2fr -Tf ~ SJ " a», 



Pu = 



fenAi^p-f^]- 



(139) 
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Some idea of the behavior of the series (134) for d n (y, a) may be 
gained from Table II. Equations (127) and (128) were programmed for 
calculation on a high speed digital computer. The table lists results for 
the typical case x = 30, a = 25, and y = -5. Here Term 2n = 
f (x)p nQ x' n , Term 2n+1 = fi(x)p ni x' n , and S m = t + h + • • • + t m . 
The "exact" value, 381.02, was calculated by using the recurrence 
relation for the Poisson-Charlier polynomials. 

No study was made to decide whether the relatively large value of 
Term 7 results from accumulated round-off error (an accuracy of 1 
part in 10 7 was used) or whether the asymptotic series actually starts 
its divergence around m = 5 or 6. 

When r is near unity, v x and t x are small and the individual terms 
in the expressions (127) and (128) for Pjo, p Jlt become large. In this 
case considerable cancellation occurs, and a high degree of precision in 
the calculations is required to obtain accurate values of Pjo and p n . 

An asymptotic series (nonuniform) which is useful when r — 1 is 
small may be obtained by a variation of the classical method which is 
sometimes used in cases of this sort. Instead of using an expansion 
about both t = and t = t x , which is done (in effect) in obtaining the 
uniform asymptotic expansion, an expansion is made only about t = 
0. Thus, the exponent xh(t) may be written as 

xh{t) = [-x(r - \)i + xt 2 /2) + (xt 2 )(t/S + i 2 /4 + • • •)• 

Changing the variable of integration from t to u = t (x/2) * and as- 
suming that r — 1 is so small that z = (r— 1) (x/2) v > is 0(1) gives 



exp [xh(t)] = exp [-2zu + w 2 ] exp [u 2 (2t/3 + 2/ 2 /4 + • • ■)] 



(140) 



= exp [-2zu + u 2 ] £ (2/.r) n/2 E b m „u 2m+n . 
The last series is obtained from equation (102) with $ = t = u (2/.t.) 1/j , 





Table II 


— Partial Sums for d 20 (— 5, 25 
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m 


Term,,, 


Sm 


m 


Ternim 


Sm 



1 
2 
o 


276.62 

S2 . 82 

20.19 

1.39 


276.62 
359.44 
379.63 
3S1.01 


4 
5 
6 

7 


0.015 

-o.oos 

0.008 
-0.144 


381.03 
381.02 
381.03 
380.88 



Exact 4- 9 (-5, 25) = 381.02 
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y = u 2 , and a„ = 2/(n+2). The coefficient 6 mn is computed from the 
a„'sand (100). 

Substituting (140) in the integral (133) for d n (y, a) leads to 



where 



d n (y, a) ~ 2 (2/aO (B+ " +1)/2 E b mn K am+n (141) 

i f e+ "° 

K n = r— . / w n f " exp (u~ — 2zu (hi) 
2iri J c -,oo 

= G(n + y, z) = G[n + y, (r - l)(.r/2) 5 J. 
The recurrence relation 

K n+1 = zK n - "^ £„-, (142) 

permits K-m + n in (141) to be expressed as a linear function of K and 
Ki. Equation (142) is obtained by integrating the derivative 

4- u n+ " exp (u 2 - 2zu) = [(n + y)u l + 2u - 2z]u n+u exp (u 2 - 2zu) 
du 

As r -» 1 the leading term, (2/xY' ,+1)/2 K , in (141) tends to the 
leading term V a {x)p 0Q in the uniform asymptotic series (134). Although 
(141) is much simpler than (134), it does not hold for nearly as wide a 
range of values of r — 1. 

APPENDIX H 

Saddle Point at Origin 

Here we are concerned with the leading term when X t± 1 and there 
are two saddle points, one at t = and the other at t = t t . We assume 
that h is real and positive, and that h(t) is real on the real axis. Further- 
more, we assume fe, < 0, K" < 0, /ij 1 " > so that the saddle point 
at t x is lower than the one at 0, and the paths of steepest descent at 
and U are parallel to the real and imaginary axes, respectively. A cut 
extends from to — <» along the negative real axis. 

The paths of integration U and L are taken to run in from co 
exp(— iir/3), cross the positive real axis in the critical region, and the 
run out to co exp(^/3). Example (Hi) of Appendix B leads us to 
choose 

F(v) = 2v 3 - 3^, vl - -I h (143) 
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with arg V\ = 0. From equation (46), the uniform asymptotic ex- 
pansion is of the form 

./ ~ V (x) E p n0 af " + VAx) Z V**" + V&) Z Pn2X~ n (144) 

,, = U n = n = 

where 

7,(3;) = J v l ^~ l exp [atfty)] dv, I = 0, 1, 2. 

Expanding exp(— SxviV 2 ) and integrating tcrmwise with the help of 

[ u p exp (w 3 ) rfw - 2ri/\ 3r (^-^)j 



gives 



V (x) = 2? (2: 



\-X/3 



E 



(-3 yi ) n (^/4)' 



n! r 1 - 



X + 2w 



(145) 



from which Vi(x) may be obtained by replacing A by A + I. When A = 
1, Vq(x) reduces to the product of an Airy function and an exponential. 
Setting n = in the integral (47) for P»(v) gives 

PoW = Poo + VoiV + po 2 y 2 

" 2t* i c nf; f(f - »,) 

- /(o) + /'(o> + [/(»,) - /(o) - ^/'(o)]^: 8 . 

The values of C> tf°» C appearing in /(0), f (0), /(i>0 are 



df (146) 



l 



= H3g 3 

L Ai a) J ' 



, 12 - KH 



(3) .(1)3 



c = 



(2) .(1) 



uHq Iq 



,u, _ [to± 1* 

1 ufj • 



(147) 



The derivatives i"', f"' are positive and nearly equal when the saddle 
points are close together. 

When n = 0, the Ursell equations (61), (62), and (63) become 



Pofyi) = aio/ftoo 

Poo = «oo//3ooo 



(148) 



«oi _ ftooi 
Poio Poio 
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The values of a ]0 , /?ioo otbained from the leading terms in the asymp- 
totic scries (54) defining the a ni 's and /?,-/„, 's give 

Po(»i) = gl n {tM x - l t?\ (149) 

Similarly, comparing the asymptotic series (64) denning the a /s and 
Point's with the scries (109) leads to expressions which give 

_ (n).(nx 
Foo — yo l o , /-. -q\ 

1_ A 

(i) ~r 



Poi = iff"' - (X+ Wo 



_3oiC G/i 




The remaining coefficient, p 2, in fo(^) may now be obtained by com- 
bining (149) and (150). The coefficients poi give the leading part of 
the desired expansion (144) for J. 
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